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Summary 

The temperatures due to frictional heating within a solid 
lubricated or coated journal bearing were analyzed by using 
a finite element method. A solid model of the shaft-bush 
tribocontact was generated with an eight-node, three- 
dimensional, first-order isoparametric heat-transfer element 
and the Patran solid modeler software. The Patmar (Patran- 
Marc) translator was used to help develop the Marc-based 
finite element program for the system; this software was used 
on the Cray X-MP super computer to perform a finite element 
analysis of the contact. The analysis was performed for various 
liner materials, for thin, hard, wear-resistant coated bearings, 
and for different geometries and thermal cooling boundary 
conditions. The analyses indicated that thermal conductivity 
of the liner or coating material is the most vital thermal 
parameter that controls the interface temperature. In addition 
to design variations, the proximity of the cooling source to 
the heat-flux-generating interface is critically important to the 
temperature control in the system. 

Introduction 

The temperature rise due to frictional heating at the interface 
of tribocontacts has attracted the attention of engineers and 
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tribologists in recent years. Friction and wear of dry contacts 
are dependent on the overall temperature increase at the 
interface. Thermomechanical effects, such as the initiation of 
thermal cracks and the low-cycle thermal-fatigue failure in dry 
bearings and seals, are dependent on the heat dissipation 
mechanism and the overall temperature rise. Floquet, Play, 
and Godet (ref. 1) used a two-dimensional Fourier transform 
method developed by Ling (ref. 2) to calculate the contact 
temperatures in dry bearings operating with a plastic liner. 
Later, Floquet and Play extended this technique to three- 
dimensional problems (ref. 3). They showed that small changes 
in the design may bring large differences in contact 
temperatures, and that to reduce the maximum temperature, 
the rubbing surfaces must be as close as possible to a controlled 
temperature heat sink and must not be shielded by any thermal 
barrier. Gecim and Winer (refs. 4 and 5) extended the Fourier 
transform method to calculate steady temperatures in a rotating 
cylinder that had multiple surface heat sources and a variety 
of geometries and thermal boundary conditions. 

In recent years, high-temperature applications brought about 
the need for hard, wear-resistant coatings (viz., A1 2 0 3 , and 
Zr0 2 ) to be used as thermal barriers in mechanical 
components. When these components are subjected to heavily 
loaded sliding conditions, “heat checking” often occurs; that 
is, cracking in the neighborhood of the contact zone due to 
the combination of thermal heating and mechanical load. This 
often occurs in such mechanical devices as seals and brakes. 
Kennedy, et al. (refs. 6 to 8) have done thermomechanical 
analyses of dry sliding contacts, using a finite element 
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technique. In particular, a recent investigation (ref. 8) 
addressed the problems of hard, wear-resistant coatings on a 
ductile metallic substrate, as applied to mechanical face seals. 
A thermo-elasto-plastic analysis of stresses and deformations 
was used. 

Ju and Chen (ref. 9) used a Fourier transform method to 
investigate the mechanical stresses, temperature field, and 
thermal stresses due to a moving frictional load. Ju and Liu 
(refs. 10 and 11) extended this technique to study the 
phenomenon of thermomechanical cracking in layered media. 
Leroy, Floquet, and Villechaise (ref. 12) also used a fast 
Fourier transform algorithm to study thermomechanical 
behavior of layered media. 

Fourier transform techniques generally provide useful 
analytical expressions, but they are limited in applications 
dealing with irregular geometries and boundaries. Insofar as 
we know, dry bearing temperatures have not been investigated 
with the finite element technique that accommodates irregular 
boundaries. Herein, results of a finite element analysis of dry 
contact temperatures due to frictional heating at the interface 
of a shaft-bush-housing system are presented. A solid model 
of the system was generated with the Patran solid modeler 
software, and the finite element analysis was done with the 
Marc finite element software. The temperature distribution in 
the system was determined with various liner materials, coating 
materials, coating thicknesses, and bearing geometries. 

This work also shows that the method can serve as a 
computer-aided design technique to model mechanical 
components (such as the cylinder-piston of an advanced low 
heat rejection (LHR) diesel engine or Stirling engine, and 
aircraft control and undercarriage bearings) for temperature 
rise and thermal distortions under specific operating 
conditions. 


Symbols 

[B] temperature gradient interpolation matrix 

[C] capacitance matrix 

c specific heat, J/kg °C 

/ friction coefficient at the contact 

h film heat transfer coefficient, W/m 2 °C 
[K] heat transfer or coefficient matrix 

k thermal conductivity, W/m °C 

M number of finite elements 

[N] temperature interpolation matrix 
n unit vector 

p contact pressure, Pa 

Q internal heat generation 

q heat flux distribution, W/m 2 or N/cm-sec 


[RJ heat load matrix 

R bi inner radius of the bearing, m 
r number of nodes in an element 

S surface area, m2 

T temperature, °C 

t time, sec 

V rubbing velocity, m/sec 

jc, y,z Cartesian coordinates 

T surface of integration, m 2 

k thermal diffusivity, m 2 /sec 

p density, kg/m 3 

Q volume of integration, m 3 

w s angular speed of the shaft, rad/sec 

Subscripts: 

c conductance 

e element, in a global sense 

h surface convection 

m maximum 

q surface heating 

s surface 

T nodal temperatures 

Superscript: 

e element, in a local sense 

Theory 

The theory discussed herein is taken from reference 13. 

Finite Element Formulation 

The steady-state or transient heat transfer in a three- 
dimensional anisotropic solid moving with velocity V is 
governed by the Fourier heat conduction or energy equation, 
which is written as 


v -k v T — pc 



F- vT 


4- 2=0 


(1) 


where Q(x,y,z t t) is the internal heat generation rate per unit 
volume. The solution for the temperature distribution must 
satisfy the boundary conditions and initial conditions on all 
portions of the surface. The inertia conditions, which specify 
the temperature at time zero, are required for the transient 
analysis. The heat conduction boundary conditions are as follows: 
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where 


where 


T=T { 

q- n = — k vT- n 
h(T s — T e ) -n = k vT-n 


( 2 ) 


^1 prescribed temperature in the specified region 

S\ of the surface 

-k vT-n prescribed heat flux normal to the surface S 2 

k v T-n convection to exchange temperature T e from the 
surface with unknown temperature T s 

The solution domain in the solid is discretized into M finite 
elements possessing r nodes. Temperature and temperature 
gradients within the element are expressed in terms of nodal 
temperatures as 


and 


T e (x,y t z,t) = [N (x,y,z,r)] (T(f)J 


dT 

— ( x,y,z,t ) 
at 


dT 

— (x.y.z.t) 
at 




(3) 


dT 


(x,y,z,t) 


where [N] is the temperature interpolation matrix, and [B] is 
the temperature gradient interpolation matrix such that, 


and 
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Finite element equations are derived from the energy equation 
(eq. 1) by using the method of weighted residuals (Galerkin’s 
method): 


[C] + [[KJ + [K„]] [71 

jdxj 


- [Rgl + [R?) + [R*l + (Rr! (5) 


[C] = 
[K c ] = 
[K a ] = 


[R Q ) = 


pc[NJ [N] cm 


[B]'[K][B]^fi 


/i(N![N]JT 


s 3 


QlNjdSI 


;rj = qMdT 


n 

[R*l = *7;sN|jr 

Ji 3 

[Rt] = - f (9-n) [NjrfT 


( 6 ) 


The coefficient matrix [C] of the time derivative of the nodal 
temperatures is the element capacitance matrix. The coefficient 
matrices [KJ and [KJ are element conductance matrices and 
relate to conductance and convection. 

The vectors (RJ, [R g], (R^j, [R h ] are heat-load vectors 
arising from specified nodal temperatures, internal heat 
generation, specified surface heating, and surface convection, 
respectively. For steady-state analysis, equation (5) reduces to 


[[K c ] + [KJ][7] = {R q } + [R,] + (R a ) (7) 


Problem Statement and Solution Methodology 

In the present work a shaft-bush/liner-housing type of 
tribosystem was analyzed. The system, shown in figure 1, 
consists of a hollow shaft (I) rotating in a bush-liner (2) fitted 
within a stationary housing (3). Alternatively, the shaft rotates 
in a housing with a thin, hard, wear-resistant ceramic coating 
at its internal surface. There is frictional rubbing at the 
interface of the liner or coating and the shaft. Frictional heating 
at the interface is conducted to the shaft and the bearing and 
housing, thereby causing a temperature rise in both. The main 
objective of the present study was to develop a finite element 
model of the previously described system by using the Patran 
solid modeler software, and thereafter to obtain a Marc finite 
element program by using a Patran-Marc (Patmar) translator. 
Following this, we analyzed the problem for temperature 
distribution within the solid, using the Marc finite element 
software, subject to the proper boundary conditions and the 
appropriate material thermal properties. The results were 
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postprocessed with the Patran postprocessor. The boundary 
conditions for the system analyzed were 

(1) A heat-transfer coefficient h equal to 58 W/m 2 °C and 
a reference temperature of 20 °C, caused by natural 
convection 

(2) An imposed temperature of 35 °C at the outer ends of 
the shaft, because of direct cooling 

(3) An imposed temperature of 35 °C at the periphery of 
the housing, because of direct cooling 

(4) An imposed temperature of 35 °C for Case I, caused 
by direct cooling at the inner radius of the shaft 

Shaft speeds of 2.6, 6.5, and 13.0 rad/sec were used to 
calculate the heat flux generated at the interface. A nodal heat 
flux q equal to the product of the contact pressure /?, the friction 
coefficient /, and the rubbing velocity V was input to the 
interface between the shaft and bush (or coated) housing. An 
harmonic contact pressure distribution extending from ±90° 
was assumed (i.e., the journal contacts of half the bearing). 

Effects of various liner materials, coating materials, coating 
thicknesses, and geometries on the temperature distribution 
and the maximum temperature at the interface were 
investigated. 

Results and Discussion 

The maximum interface temperatures and the temperature 
distribution in the shaft-bearing system are presented herein. 
This system was composed of either a bearing with liner 
material in a housing, or a bearing with a thin, hard, wear- 
resistant coating on its inner surface. The temperature rise was 
due to frictional heating at the interface. The hollow shaft had 
different dimensions (see table I) for the three cases studied: 
(1) bearings with soft liner materials; (2) bearings with thin, 
hard coatings on the internal surface; and (3) variations in 
geometry (table II) and the cooling condition. Nodal heat flux 
generated in the contact was estimated to be equal to the 
product of the contact pressure, rubbing velocity, and friction 
coefficient. The coefficients of friction for the liner and coating 
materials ranged from 0.15 to 0.3. However, an intermediate 
value of 0.2 was assumed. Our objective was to evaluate the 
temperature rise in the contact for the same amount of heat 
flux input. A harmonic pressure distribution extending from 
±90° was assumed. 

Thus, heat flux distribution 

q = q m cos <t> 

where 

Rm “ jVPm 
and 

/ friction coefficient at the contact, assumed to be 0.2 
V rubbing velocity; R bl w s 
p m maximum contact pressure 


Case I— Bearings with Soft Liner Materials 

Thermal properties of various liner materials are shown in 
table III. The shaft and the housing were made of steel. For 
a 6.0-mm thick liner, the maximum interface temperatures 
obtained for various liner materials are shown in table IV. The 
lowest of the maximum temperatures was 108.5 °C for the 
graphitized metal liner, which had the highest thermal 
conductivity and diffusivity values. As the material’s 
conductivity and diffusivity decreased, the temperature at the 
interface increased. In the absence of convection, the thermal 
conductivity of the liner was the material’s single most 
important thermal property, influencing both the temperature 
rise at the interface and the temperature distribution as well 
(fig. 2). Table V shows the variation of maximum interface 
temperature with five different heat distribution functions for 
a resin liner material (ref. 1), while keeping all other conditions 
constant. We observed that the maximum interface temperature 
and the temperature distribution can differ significantly with 
the nature of the heat distribution in the contact region (fig. 3). 
The nature of the heat distribution depends on the contact 
pressure in the contact region, which in turn depends on the 
geometry and the type of loading. The nature of the loading 
may vary with the operating conditions. 

The effects of the magnitude of the load and the thickness 
of the liner material are shown in table VI. Not surprisingly, 
a thicker liner material can result in a lower interface 
temperature, provided the conductivity of the liner material 
is higher than the conductivity of the bush material. Thus, for 
higher loads or pV factors, a thicker graphitized metal liner 
would help reduce the temperatures. 

Case II— Bearings with a Thin, Hard Coating on the 
Internal Surface 

Thermal properties of thin, hard, coating materials are given 
in table VII. For a 1.0-mm-thick coating, the variations in 
maximum interface temperatures are shown in table VIII. As 
is evident from table VII for the thermal properties and table 
VIII for the maximum temperatures, the thermal conductivity 
of the coating material determines the magnitude of the 
maximum interface temperature. When a coating material with 
low thermal conductivity was used as a thermal barrier to 
shield the bearing, high temperatures occurred at the shaft- 
bearing interface, thereby dissipating more heat to the shaft 
and causing higher shaft temperatures. 

In order to study the effect of coating thickness, calculations 
were made for titanium nitride coatings ranging from 
0.125 to 1.0 mm thick. The resulting maximum interface 
temperatures are shown in table IX. Note that maximum 
interface temperatures decreased with decreasing coating 
thickness. This means that thicker coatings with lower 
conductivity than the substrate (such as titanium nitride on 
steel) are better thermal barriers than very thin coatings. 
Temperature distributions in the system, which result from 
using three different thicknesses of titanium nitride coating, 
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are shown in figures 4 to 6. Similar calculations for various 
thicknesses of tungsten carbide coatings showed almost no 
change in maximum interface temperature, since tungsten 
carbide has a higher thermal conductivity than the steel 
substrate. 

Case III— Variations in Geometry and the Cooling 
Condition 

Tables X and XI show the variation of maximum interface 
temperatures for a number of bearing outer radii and bearing 
lengths. Other dimensions remained the same. Higher interface 
temperatures occurred with larger bearing diameters, since 
the cooling source was located at the outer bearing diameter. 
An increase in bearing length did not alter the maximum 
temperature. For a hollow shaft that was cooled internally, 
the interface temperature was drastically reduced— from 203 
to 109 °C. The thermal properties of the materials are 
considered to be temperature-independent here. Figures 7 to 
10 show the changes in the temperature distribution in the 
system because of geometry variation and inner surface cooling 
of the hollow shaft. 

An increase in rubbing velocity or pV factor can cause a 
tremendous rise in the maximum interface temperatures. Table 
XII shows that temperatures as high as 928.0 °C can be 
obtained at pV — 5362.5 N/cm-sec with a 1.0-mm thick 
titanium nitride coating. Figure 11 illustrates the thermal 
distribution for this case. 


Summary of Results 

A Patran solid modeling package was used to generate a 
finite element model for the thermal analysis of a shaft-bush 
tribosystem. The Marc finite element software was used to 
determine the temperature distribution in the tribosystem. 
These software permit finite element modeling of bearings and 
seals so that temperature distributions can be calculated in 
specific applications. Mechanical components (viz., the 
cylinder-piston of an advanced LHR diesel engine, the Stirling 
cycle engine for space power applications, aircraft control and 
undercarriage bearings, and so on) can be modeled for 
temperature-rise and thermal distortions under specific 
operating conditions. The results of our analysis can be 
summarized as follows: 

1. If a constant friction coefficient is assumed, thermal 
conductivity of the liner or coating material is the most 
important thermal property for controlling temperatures in 
tribocontacts. 

2. Design changes in bearing diameter and length can result 
in significant differences in the interface temperature. 

3. Interface temperatures can be substantially reduced by 
cooling the inner surface of the hollow shaft. 

4. Designs that can provide a cooling medium in close 
proximity to the heat-flux-generating interface would be most 
effective in reducing the temperature rise in the contact. 
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TABLE I.— SHAFT SPECIFICATION AND BEARING 
DIMENSIONS FOR THREE CASES 


Case 

Shaft 

Bearing 


Outer 

Inner 

Length, 

Outer 

Inner 

Length, 


radius, 

radius, 

cm 

radius, 

radius. 

cm 


cm 

cm 


cm 

cm 


a I and a II 

2.5 

1.5 

8.4 

6.0 

2.5 

1.2 

•>111 

2.5 

1.5 

16.0 

5.0 

2.5 

1.6 


a Bearing is set off shaft to right by 1 .2 cm. 
^Bearing is centrally located on shaft. 


TABLE II.— GEOMETRY 
VARIATIONS FOR 
CASE III 


TABLE IV. — MAXIMUM 
INTERFACE TEMPERATURES 
FOR VARIOUS LINER 
MATERIALS 


[Heat flux distribution q = <? m cos0; 
pV = 1072.5 N/cm-sec; liner 
thickness = 6.0 mm.] 


Liner material 

Maximum 

temperature, 

°C 

Graphitized metal 

108.5 

PTFE-bronze 

118.6 

Carbon graphite 

147.0 

PTFE-thermoset 

152.8 

PTFE-ceramic 

152.7 

Mo 2 S-nylon 

152.8 

Resin 

152.8 


Bearing 

Bearing outer 

length, 

radius, 

cm 

cm 

1.6 

3.2 

1.6 

5.0 

1.6 

8.0 

3.2 

5.0 


TABLE V.-MAXIMUM INTERFACE 
TEMPERATURES FOR VARIOUS HEAT 
FLUX DISTRIBUTION FUNCTIONS 


[Resin liner thickness = 6.0 mm; pV = 1072.5 
N/cm-sec.] 


Heat flux distribution, 

q, 

N/cm-sec 

Maximum 

temperature, 

°C 

Maximum 

temperature 

location, 

deg 

4 m COS <f) 

152.8 

0 

<7 m *(4/7T)COS 2 0 

179.7 

0 

^ m *(2/7r)(l +sin</>) 

169.7 

±90 

^•(2/tt) 

114.0 

0 

q m * (2/ 7r)( 1 sin0) 

99.6 

0 


TABLE IIL— THERMAL PROPERTIES OF LINER MATERIALS 


[From ref. 14.] 


Liner Material 

Specific 

heat, 

c, 

J/kg °C 

Thermal 

conductivity, 

K 

W/m °C 

Density, 

P, 

kg/m 3 

Diffusivity, 

KlpC , 

m 2 /sec 

Graphitized metal 

377.0 

127.4 

7039.0 

4.8 X 10-5 

PTFE-bronze 

394.0 

42.5 

6501.6 

1.66 x 10-5 

Carbon graphite 

754.0 

3.35 

2101.7 

2.1 x 10-5 

PTFE-thermoset 

942.8 

.263 

2293.2 

1.1 X 10" 6 

PTFE-ceramic 

838.0 

.337 

2376.0 

1.69 X 10- 7 

Mo 2 S-nylon 

188.6 

.25 

1198.0 

1.22 x 10" 7 

Resin 

1214.0 

.233 

1200.0 

1.6 x 10" 7 


TABLE VI.— MAXIMUM INTERFACE TEMPERATURES 
FOR pV FACTORS AND LINER THICKNESSES 

[Heat flux distribution q = q m cos </>.] 

(b) Graphitized metal liner; 

(a) Resin liner, 6.0 mm thick pV = 1072.5 N/cm-sec 


Liner 

Temperature, 

thickness, 

°C 

mm 


6.0 

108.5 

12.0 

103.5 

18.0 

99.8 


p m V factor, 

Temperature, 

N/cm-sec 

°C 

1072.50 

152.8 

536.25 

93.9 

107.25 

46.8 
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TABLE VII. -THERMAL PROPERTIES OF COATING 
MATERIALS AND SUBSTRATE 


TABLE X.— MAXIMUM INTERFACE 
TEMPERATURES FOR VARIOUS 
BEARING DIMENSIONS 


Material 

Specific 

heat, 

c, 

J/kg °C 

Thermal 

conductivity, 

K 

W/m c C 

Density, 

P, 

kg/m 3 

Diffusivity 
tdpc , 
m 2 /sec 

Aluminum 

869.8 

267.7 

2 

699.0 

1.136x10-5 

Zirconium 

17 769.0 

131.9 

6 

513.0 

1.142x10-6 

Silicon carbide 

711.0 

104.4 

2 

997.0 

4.9x10-5 

Tungsten carbide 

100.0 

100.0 

15 

800.0 

5.0x10-5 

Aluminum oxide 

1 039.0 

20.73 

3 

997.0 

4.99x10-5 

Stellite 

421.9 

9.7 

8 

300.0 

2.77x10-5 

Titanium nitride 

590.0 

4.0 

5 

200.0 

1.3x10-5 

Zirconia 3 

400.0 

1.8 

5 

750.0 

7.8x10-5 

Steel 

460.0 

46.0 

7 

800.0 

1.28x10-5 


a Partially stabilized. 


TABLE VIII. — MAXIMUM 
INTERFACE TEMPERATURES 
FOR COATING MATERIALS 


[Coating thickness = 1.0 mm; heat 
flux distribution q = q m cos<f>; 
pV = 1072.5 N/cm-sec.] 


Coating material 

Maximum 

temperature, 

°C 

Aluminum 

176.8 

Zirconium 

181.1 

Silicon carbide 

182.2 

Tungsten carbide 

182.4 

Aluminum oxide 

191.3 

Stellite 

199.5 

Titanium nitride 

213.0 

Zirconia 3 

234.7 


a Partially stabilized 


TABLE IX.— MAXIMUM 
INTERFACE 
TEMPERATURES FOR 
VARIOUS TITANIUM 
NITRIDE COATING 
THICKNESSES 


[Heat flux distribution 
q = q m cos<f>; pV 
= 1072.5 N/cm-sec.] 


Coating 

thickness, 

mm 

Maximum 

temperature, 

°C 

1.0 

213.0 

.75 

208.0 

.5 

201.0 

.25 

194.0 

.125 

190.0 


[Coating, 1-mm-thick TiN; heat 
flux distribution q = q m cos<t>\ 
pV = 1072.5 N/cm-sec.] 


Bearing and housing 

Maximum 

outer radius, 

temperature, 

cm 

°C 

3.8 

182.0 

5.0 

202.1 

8.0 

228.1 


TABLE XI.— MAXIMUM 
INTERFACE TEMPERATURES 
FOR VARIOUS BEARING 
LENGTHS AND COOLING 
CONDITIONS 


[Coating, 1-mm-thick TiN; 
pV = 1072.5 N/cm-sec; bearing 
outer radius, 5 cm.] 


Bearing length, 

Maximum 

cm 

temperature, 


°Q 

1.6 

202.1 

3.2 

202.5 

a 3.2 

109.0 


a With shaft inner surface cooling. 


TABLE XII.— MAXIMUM 
INTERFACE 
TEMPERATURES FOR 
VARIOUS pV FACTORS 


[Coating, 1-mm-thick TiN; 
heat flux distribution 
q = q m cos</>.] 


pV factor, 

Maximum 

N/cm-sec 

temperature, 


°C 

1072.5 

213.0 

2681.25 

481.7 

5362.5 

928.0 
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(a) Section view. 

(b) Finite element model. 

Figure 1.— Shaft-bush bearing tribosystem. 
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Figure 2.— Tribosystem temperature contours for case I. Bearing has 6.0-mm-thick soft resin liner. 
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(a) Heat distribution function q = q m { 2/n) (1 + sin <£). 

Figure 3.— Tribosystem temperature contours for case I. Bearing has 6-mm-thick soft resin liner and various heat distribution functions. 
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(b) Heat distribution function q = q m ( 2 /tt). 
Figure 3.— Continued. 
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(c) Heat distribution function q = q m ( 2/n) (l-sin$). 
Figure 3.— Concluded. 
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Figure 4.— T ribosystem temperature contours for case II. Bearing has 1.0-mm-thick TiN coating. 
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Figure 5. — Tribosystem temperature contours for case II. Bearing has 0.5-mm-thick TiN coating. 




Figure 6.— Tribosystem temperature contours for case II. Bearing has 0. 125-mm-thick TiN coating. 


Temperature, 
°C 

190 

179 

169 

159 

148 

138 

128 

117 

107 

96.8 

86.5 

76.2 

65.9 

55.6 

45.3 
35.0 



Temperature, 

°C 



Figure 7.— Tribosystem temperature contours for case III. Bearing has 1.0-mm-thick TiN coating; bearing outer radius = 8.0 cm. 
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Figure 8.-Tribosystem temperature contours for case III. Bearing has 1.0-mm-thick TiN coating; bearing outer radius = 3.8 cm 
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Figure 9.— Tribosystem temperature contours for case III. Bearing has 1.0-mm-thick TiN coating; bearing outer radius 
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Figure 10.— Tribosystem temperature contours for case III. Bearing has 1.0-mm-thick TiN coating and shaft inner surface cooling; bearing outer 
radius = 5 cm; bearing length = 3.2 cm. 
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Figure 1 1.— Tribosystem temperature contours for case II. Bearing has 1.0-mm-thick TiN coating; pV factor = 5362.5 N/cm-sec. 
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